syms x
int((x^4*(1-x)^4)/(1+x^2),0,1)
f=@(x)(x.^4.*(1-x).^4)./(1+x.^2);
format long
quad(f,0,1)
quadl(f,0,1)
format default